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ABSTRACT 

Context. The space experiment CoRoT {Convection, Rotation and Transits) has recently detected a transiting hot Jupiter in orbit 
around a moderately active F-type main-sequence star (CoRoT-Exo-4a). This planetary system is of particular interest because it has 
an orbital period of 9.202 days, the second longest one among the transiting planets known to date. 

Aims. We study the surface rotation and the activity of the host star during an uninterrupted sequence of optical observations of 58 
days. 

Methods. Our approach is based on a maximum entropy spot modelling technique extensively tested by modelling the variation of the 
total solar irradiance. Recently, it has been successfully applied to model the light curve of another active star with a transiting planet 
observed by CoRoT, i.e., CoRoT-Exo-2a. It assumes that stellar active regions consist of cool spots and bright faculae, analogous to 
sunspots and solar photospheric faculae, whose visibility is modulated by stellar rotation. 

Results. The modelling of the light curve of CoRoT-Exo-4a reveals three main active longitudes with lifetimes between ~ 30 and 
~ 60 days that rotate quasi-synchronously with the orbital motion of the planet. The different rotation rates of the active longitudes 
are interpreted in terms of surface dilferential rotation and a lower limit of 0.057 ± 0.015 is derived for its relative amplitude. The 
enhancement of activity observed close to the subplanetary longitude suggests a magnetic star-planet interaction, although the short 
duration of the time series prevents us from drawing definite conclusions. 

Conclusions. The present work confirms the quasi-synchronicity between stellar rotation and planetary orbital motion in the CoRoT- 
Exo-4 system and provides for the first time a lower limit for the surface differential rotation of the star. This information can be 
important in trying to understand the formation and evolution of this highly interesting planetary system. Moreover, there is an 
indication for a possible star-planet magnetic interaction that needs to be confirmed by future studies. 

Key words, stars; magnetic fields - stars: late-type ~ stars: activity - stars: rotation - planetary systems - stars: individual (CoRoT- 
Exo-4a) 



1. Introduction 

CoRoT is a photometric space experiment devoted to asteroseis- 
mology and the search for extrasolar planets by the method of 
transits teaghn et al.l l2006V It has recently discovered CoRoT- 
Exo-4b, a Jupiter-sized planet transiting across the disc of an 
F-type main-sequen ce star with an orbital period of 9.202 days 
jAigrain et alJl200 8V This is the second longest period among 
the transiting planetary systems known to date, putting CoRoT- 
Exo-4b in a region of the mass-period para meter space that was 
previously empty (cf. iMoutou et al.ll2008b . Moreover, the out- 
of-transit light curve shows a modulation with an amplitude of a 
few 0.001 mag that can be attributed to photospheric brightness 



* Based on observations obtained with CoRoT, a space project op- 
erated by the French Space Agency, CNES, with partecipation of the 
Science Programme of ESA, ESTEC/RSSD, Austria, Belgium, Brazil, 
Germany, and Spain. 



inhomogeneities carried into and out of view by the rotation of 
the star. Given its late spectral type, those inhomogeneities can 
be considered analogous to cool spots and bright faculae ob- 
served in the Sun, owing their existence to photospheric mag- 
netic fields. 

The active regions of CoRoT-Exo-4a are sufficiently stable 
to allow an estimate of its rotation period through an autocorre- 
lation analysis, yielding a period of 8.87 + 1.12 days. This re- 
sult indicates that the stellar rotation and the orbital motion of 
the planet are quasi-synchronized, which cannot be explained 
on the basis of the tidal theory given the large separation be- 
tween the star and the planet (a/Rs mr = 17.36 ± 0.25) and its 
low mass (0.72 + 0.08 Mjup, cf., lAigrain et al]|2008l) . Other 
synchronous or quasi-synchronous star-planet systems have re- 
cently been foun d, consisting of quite massive planets orbit - 
ing F-type stars (ICatala etal.ll2007t iMcCullough etaH l2009l) . 
Among those systems, r Bootis has been recently studied in 
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some detail revealing a surface differential rotation with a rela- 
tive amplitude AQ/Q =s 0.15 -0.18 and hints for an activity cycle 
of a few years dCatala et al J'20Q7; Donati et al. 2008). However, 
those systems host more massive and closer planets that may ex- 
plain the synchronization on the basis of tidal effects, assuming 
that only the outer stellar convection zone, having a mass of only 
a few 0.001 of the stellar mass, is synchronized (cf. iDonati et alj 
I2OO8). On the other hand, the estimated time scale for the syn- 
chronization of the outer convective envelope of CoRoT-Exo-4a 
is =s 350 Gyr, assuming a mass of 0.001 of the stellar mass for 
the envelope. 

In the present work, we present a modelling of the out-of- 
transit light curve of CoRoT-Exo-4a following the approach al- 
ready applied for another star hosting a transiting hot Jupiter, i.e., 
CoRoT-Exo-2a. It is based on a solar analogy, i.e., it assumes 
that the stellar active regions responsibl e for the flux modu lation 
consist of cool spots and bright faculae dLanza et~al ] l2009h . 

Late-type stars accompanied by hot Jupiters sometimes show 
active regions rotating with the orbital period of the planet in- 
stead of the stellar rotation period. The phenomenon has been 
detected in some seasons in v Andromedae and HD 179949 that 
have shown chromospheric ho t spots leading the planet by ~ 
170° and ~ 70°, respectively dShkolnik et al.| [2008). Moreover, 
the synchronous system r Boo has shown some evidence for an 
active region rotating with the orbital period of it s planet but 
leadin g it by ^ 70° during the period 2001-2005 (Wal ker etal] 
I2OO8I) . These phenomena suggest a magnetic interacti on be - 
tween the star and its hot Jupiter, as discussed bv lLanzal (l2008h . 
In the case of CoRoT-Exo-4a, the short duration of the obser- 
vations prevents us from deriving definite conclusions about a 
possible star-planet interaction, but the longitude distributions 
of the active regions may be compared with the orbital position 
of the planet to see whether there is some analogy, say, with t 
Boo. 

A detailed study of stellar activity by means of Doppler 
Imaging techniques is not feasible owing to the small projected 
rotation velocity of CoRoT-Exo-4a (vsin; = 6.4 + 1.0 km s"') 
and its faintness (R ~ 13.5). Therefore, a spot modelling ap- 
proach is the most suitable method to investigate the magnetic 
activity and rotation of this very interesting object. 



The light curve extracted from the CoRoT Mission Archiv43 
was further corrected for the effect of a hot pixel in the blue 
channel, as explained by Aigrain et al. (2008). To increase the 
signal-to-noise ratio and reduce residual systematic effects pos- 
sibly present in individual channels, we summed up the flux in 
the red, green and blue channels to obtain a light curve in a spec- 
tral range extending from 300 to 1 100 nm. The oversampled sec- 
tion of this white-band light curve was rebinned with a regular 
512 s sampling, to obtain an evenly sampled time series. It was 
further cleaned by applying a moving median 5-cr clipping fil- 
ter that allowe d us to identify an d discard residual outliers, as 
described by Aigr ain et al.l (|2008). The standard deviation of the 
observations, averaged over 512 s bins, is 8.9 x 10""* in rela- 
tive flux units. The contamination of the photometric aperture 
by stars other than CoRoT-Exo-4 was only 0.3 + 0.1 percent of 
the median of the measured flux and it was subtracted to avoid 
dilution of the stellar light variations. Tran sits were removed b y 
means of the ephemeris and parameters of iAigrain et alj (l2008l) . 

The flux variations related to stellar activity have time scales 
of the order of a day or longer Therefore, we rebinned the out- 
of-transit light curve by computing normal points along each or- 
bital period of the satellite (6184 s). This has the advantage of 
removing tiny systematic variations that ma y still be associated 
with the orbital motion of the satellite (cf. lAlonso et al.l 120081 : 
lAuvergne et ani20b 9). Since the orbital period is not a multiple 
of the individual exposures of 32 or 512 s, we performed a lin- 
ear interpolation of the flux variation along the last 40 s of each 
orbital period to compute the normal points. In such a way we 
obtained a light curve consisting of 684 normal points ranging 
from HJD 2454135.0918 to HJD 2454192.7744, i.e., with a du- 
ration of 57.6826 days. 

The hg ht curve shows a l ong-term decreasing trend (cf. 
Fig. 2 in 'Aigrain et al.' '2008) that may be of instrumental 
origin ( Auvergne et al. 2009). Since it would produce only a 
monotonous increase of the spotted area in our model along the 
~ 60 days of the present observations, we decided to subtract it 
before modelling the light curve by fitting a straight line to the 
data. Finally, the de-trended light curve was normalized at its 
maximum flux value observed at HJD - 2455737.2730, that we 
assumed to represent the unspotted flux level of the star, whose 
true value is unknown. 



2. Observations 

CoRoT-Exo-4a was observed during the CoRoT Initial Run for 
58 days, starting from 6 February 2007. The time sampling was 
initially of 5 12 s, then reduced to 32 s during the last three weeks 
after the detection of the transits by the so-called CoRoT alarm- 
mode. CoRoT performed aperture photometry with a fixed mask. 
Since the star was brighter than R - 14.5, the flux was split along 
detector column boundaries into broad-band red, green and blue 
channels. 

The obse r vation s and data processing are described by 
lAigrain et alj ( 1200 8') to whom we refer the reader for de- 
tails. The reduction pipeline applied corrections for the back- 
ground and the pointing jitter of the satellite which was par- 
ticularly relevant during ingress and engress from the Earth 
shadow. Measurements during the crossing of the South Atlantic 
Anomaly of the Earth magnetic field, which amounted to about 
15-20 percent of a satellite orbit, were discarded. More infor- 
mation on the instrumen t , its op eration and performance can be 
found in lAuvergne et al.l (l2009l) . 



3. Spot modelling 

The reconstruction of the surface brightness distribution from 
the rotational modulation of the stellar flux is an ill-posed prob- 
lem, because the variation of the flux vs. rotational phase con- 
tains information only on the distribution of the brightness inho- 
mogeneities vs. longitude. The integration over the stellar disc 
cancels any latitudinal information, particularly when the incli- 
nation of the rotation axis along the line of sight is 90°, as in 
the present case (see Sect.|4]and iLanza et"ani2009[) . Therefore, 
we need to include a priori information in the light curve in- 
version process to obtain a unique and stable map. This is done 
by computing a Maximum Entropy (hereinafter ME) map which 
has been proved to give the best representation in the case of the 
Sun (cf. iLanza et alj|2007l) . As a matter of fact, a sequence of 
seven maps is obtained for CoRoT-Exo-4a, covering its succes- 
sive rotations along the 58 days of the observations. They will 
allow us to study the evolution of the active longitudes over the 
surface of the star Readers not interested in the details of our 
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spot modelling approach may skip the remainder of this Section 
and go on to Sect. |4] (or directly to Sect.|5]) • 

In our model the star is subdivided into a large number of 
surface elements, in our case 200 squares of side 18°, with each 
element containing unperturbed photosphere, cool spots and fac- 
ular areas. The fraction of an element covered by cool spots is 
indicated by the filling factor /, the fractional area of the faculae 
is Qf and the fractional area of the unperturbed photosphere is 
1 - (2 + !)/■ We fit the light curve by changing the value of / 
over the surface of the star, while Q is held uniform and constant. 
Even fixing the rotatio n period, the inc lination, and the spot and 
facular contrasts (see Lanza et aLll2007l for details), the model 
has 200 free parameters and suffers from non-uniqueness and 
instability. To find a unique and stable spot map, we apply maxi - 
mum entropy regularization, as described in lLanza et alj (|2007[) . 
by minimizing a functional which is a linear combination of 
the;\f^ and the entropy functional 5 , i.e.: 

&^x\f)-^s(f\ (1) 

where / is the vector of the filling factors of the surface ele- 
ments, /t > a Lagrangian multiplier determining the trade-off 
between light c urve fitting and regu larization, and the expression 
of S is given in lLanza et al.l ( Il998h . The entropy functional S is 
constructed in such a way that it attains its maximum value when 
the star is immaculate. Therefore, by increasing the Lagrangian 
multiplier A, we increase the weight of S in the model and the 
area of the spots is progressively reduced. This gives rise to sys- 
tematically negative residuals between the observations and the 
best fit model when A > Q. The optimal value of A depends on 
the information content of the light curve, that in turn depends 
on the ratio of the amplitude of its rotational modulation to the 
average standard deviation of its normal points. To optimize the 
extraction of this information, we generalize the cri terium to fix 
the La grangian multiplier previously introduced in iLanza et"aLl 
(l2009h . For a given value of A, we compute the mean of the 
residuals for the regularized best fit yUreg and compare it with 

eo = (To/ V^, i.e., the standard error of the residuals, where ctq 
is their standard deviation, obtained in the case of the unregu- 
larized best fit (i.e., for A - 0), and is the number of normal 
points in each fitted subset of the light curve of duration Aff (see 
below). We iterate on the value of A until l/iregl = P^o, where (3 
is a numerical factor that will be fixed a posteriori according 
to the requisites of an acceptable fit and a regular overall evo- 
lution of the spots (cf. Sect. |5]l. In the case of the light curve 
of CoRoT-Exo-2a, the rotational modulation has an amplitude 
of ~ 0.06 mag, which corresponds to a signal-to-noise ratio of 
~ 300 for a standard deviation of a; 2 x 10"'* mag. In the present 
case, the amplitude of the rotational modulation of CoRoT-Exo- 
4a is only ~ 0.006 mag and the standard deviation of the points 
is somewhat greater because the star is fainter, giving a signal- 
to-noise ratio of ~ 20. Therefore, while /3 - I was found to be 
adequate for CoRoT-Exo-2a, /? > 1 is required to appropriately 
reconstruct ME maps of CoRoT-Exo-4a. 

In the case of the Sun, assuming a fixed distribution of the 
filling factor, it is possible to obtain a good fit of the irradiance 
changes only for a limited time interval Aff, not exceeding 14 
days that is the lifetime of the largest sunspot groups dominat- 
ing the irradiance variation. In the case of other active stars, the 
value of Aff must be determined from the observations them- 
selves, looking for the maximum data extension that allows for a 
good fit with the applied model (see Sect.|4]for CoRoT-Exo-4a). 

The optimal values of the spot and facular contrasts, and of 
the facular-to-spotted area ratio Q in stellar active regions are 



unknown a priori. In our model the facular contrast Cf and the 
parameter Q enter as the product CfQ, hence we can fix Cf and 
vary Q, estimating its best value by minimizing the of the 
model, as will be shown in Sect. H] Since the number of free 
parameters of the ME model is lar ge, for this specific a pplication 
we shall make use of the model of iLanza et a n (l2003h which fits 
the light curve by assuming only three active regions to model 
the rotational modulation of the flux plus a uniformly distributed 
background to account for the variations of the mean light level. 
This procedure is the same adopted to fix the value of Q in the 
case of CoRoT-Exo-2a (cf. iLanza et"al ] l2009t) . 

As in the case of CoRoT-Exo-2a, we shall assume an inclina- 
tion of the rotation axis of CoRoT-Exo-4a of 90°, which implies 
that our model cannot provide information on the latitude dis- 
tribution of the active regions. Such a limitation cannot be over- 
come by assuming an inclination / < 90° in the model, owing to 
the relatively low information content of the present light curve. 
Specifically, when we assume / < 90°, we find that the ME regu- 
larization virtually puts all the spots at the sub-observer latitude 
(i.e., 90° - /) to minimize their area and maximize the entropy. 
Therefore, we are limited to map only the distribution of the ac- 
tive regions vs. longitu de, that can b e done with a resolution of 
at least ~ 50° (cf. .Lanza et alj|2'007l [2009) . Note that our igno- 
rance of the true value of the facular contribution to the light 
variations may produce systematic errors in the active region 
l ongitu des derived by our model, as discussed by .Lanza et all 
( |2007|) in the case of the Sun. 

4. Model parameters 

The fu nda mental stellar paramete rs are taken froml Aigrain et al.l 
( I2008h and iMoutou et all (l2008h and are Usted in Table [1] A 
quadratic limb-darkening law is adopted for the stellar photo- 
sphere, viz., cc (flp + bpfi + Cp/v^), where /(//) is the spe- 
cific intensity at the limb position ju s cos 9, with 9 being the 
angle bet ween the normal to a su rface element and the line of 
sight (cf. iLanza et"alll2003L l2007h. The hm b-darkening param- 
eters have been derived from'Ku ruczl (|200C|) model atmospheres 
for Tgff = 6190 K, logg = 4.41 (cm s"^) and solar abun- 
dances, adoptirig the C oRoT white-band transmission profile in 
lAuvergne et anl2009l) . 

The polar flattening of the star due to the centrifugal potential 
is computed in the Roche approximation with a rotation period 
of 9.202 days (see below). The relative difference between the 
equatorial and the polar radii is e = 1.10 x lO"'* which may 
induce a flux variation about two orders of magnitude smaller for 
a spot coverage of ~ 1 percent as a consequence of the gravity 
darkening of the equatorial regions of the star. 

The inclination of the stellar rotation axis is difficult to con- 
strain through the observation of the Rossiter-McLaughlin ef- 
fect because of the small v sin / of the star and its intrinsic line 
profile variations due to stellar activity (Moutouetal. 20081). 
Nevertheless, we shall assume that the stellar rotation axis is 
normal to the orbital plane of t he planet, i.e., at an inclination 
of 90° from the line of sight (cf. lAigrain etal]|2 008). 

The rotation period adopted for our spot modelling is equal 
to the orbital period of the planet. This will allow us to check the 
synchronization of stellar rotation to the planetary orbit through 
the study of the longitude drift of the active regions versus time 
(cf. Sect. |5]l. Another advantage of this choice is that the sub- 
planetary longitude is fixed in our reference frame, which al- 
lows us to investigate straightforwardly the possible associa- 
tion between active longitudes and orbital position of the planet. 
Nevertheless, considering that the autocorrelation analysis of the 
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Fig. 1. The ratio of the of the composite best fit to the en- 
tire time series of CoRoT-Exo-4a to its minimum value vs. the 
parameter Q, i.e., the ratio of the area of the faculae to that of 
the cool spots in active regions. The horizontal dashed line in- 
dicates the 99 percent confidence level ior x^/x^in^ determining 
the interval of the acceptable Q values. 

Table 1. Parameters adopted for the spot modelhng of CoRoT- 
Exo-4a. 



Parameter 




Ref." 


Star Mass (Mq) 


1.16 


AOS 


Star Radius (Ro) 


1.17 


AOS 


Ta (K) 


6190 


M08 




0.334 


L09 




1.032 


L09 


Cp 


-0.381 


L09 


/'ro, (d) 


9.20205 


AOS 


6 


1.10 X 10 


L09 


Inclination (deg) 


90.0 


AOS 




0.681 


L09 


Cf 


0.115 


L04 


Q 


4.5 


LOS 


Aff (d) 


8.24037 


L09 



" References: AOS: Aigrain al. (2008); M08: Moutou et al. (2008); L04: 
Lanza et al. (2004); L09: present study. 



light mod ulation gives a rotation period of 8.87 + 1.12 days 
jAigrain et al . 2008), it is worth investigating also spot models 
with a rotation period of 8.87 days. 

The maximum time interval that our model can accurately fit 
with a fixed distribution of active regions, has been determined 
by means of a series of tests and has been found to be of ~ 8.2 
days. Therefore, we subdivided the time series into seven equal 
segments of duration Aff = 8.24037 days, each of which was 
modelled with a fixed active region pattern. 

To compute the spot contrast, we adopted the same mean 
temperature difference derived in the ca se of sunspot group s 
from their bolometric contrast, i.e. 560 K ( Chapman et al.ll994h . 
In other words, we assumed a spot effective temperature of 5630 
K, y ielding a constrast Cj - 0.68 1 in the CoRoT passband (cf. 
iLan za et al. 2007). A different spot contrast changes the absolute 
spot coverages, but does not significantly affect thei r longitudes 
and their time evolution, as discussed in detail bv iLanza et al.l 



(l2009h . Th e facular contrast w as assumed to be solar-like with 
Cf = 0.115 ('Lanza et al.ll2004 . 

The best value of the area ratio Q between the faculae and the 
spots in e ach active region h as been estimated by means of the 
model bv lLanza etal] (I2003L cf. Sect. In Fig. [1] we plot the 
ratio of the total of the composite best fit to the entire 

time series to its minimum value . , versus Q. The horizontal 

^ mm' *^ 

dashed line indicates the 9 9 percent confidence level as derived 
from the F-statistics (e.g., iLampton et al.|[T976l) . The best value 
of Q turns out to be 2 = 4.5, with an acceptable range extending 
from ~ to ~ 6.5. For comparison, the best value in the case of 
the Sun is Qq = 9, indicating a lower relative contribution of the 
faculae to the total light variation in CoRoT-Exo-4a. The ampli- 
tude of the rotational modulation of our star is ~ 0.006 mag, i.e., 
2-3 times that of the Sun at the maximum of the eleven-year 
cycle (see Sect. |5j. This suggests that CoRoT-Exo-4a is more 
active than the Sun, which may account for the red uced facular 
contrib utio n to its hght variations, as suggested bv iRadick et al.l 
(119981) and lLockwood et all (l2007h . 

5. Results 

The composite best fit obtained with the ME regularization is 
shown in Fig. |2] together with the residuals. The systematic 
negative values of the residuals are a consequence of the ME 
regularization that tends to reduce the spotted area as much as 
possible, thus systematically driving the fit above the observa- 
tions (cf. Sect.|3]l. 

The residuals of the unregularized composite best fit have a 
mean of 7.845 x 10"^ and a standard deviation ctq = 2.832 x 10""^ 
in relative flux units, i.e., assuming the maximum observed flux 
as the unit of measure. Note that in the absence of regulariza- 
tion the standard deviation of the residuals is larger t han ex- 
pecte d from the photon shot noise, i.e., 1 .61x10 '* (Aigrai iTet al.l 
l2008l) . This may arise from intrinsic stellar microvariability on 
short-term time scales (up to a few days) or may be due to some 
residual instrumental effects. On the average, we have ~ 100 
normal points per fitted subset. The Lagrangian multipliers for 
the regularized ME models have been fixed to obtain a mean 
of the residuals of the composite fit /ireg = -2.731 x 10""^, cor- 
responding to lyUregl - 9.6 fQ. The optimal value /? - 9.6 was 
established by considering the information content of the maps 
obtained for different values of /3 and minimizing it (i.e., maxi- 
mizing the entropy) as much as possible, while still retaining an 
average acceptable fit. Smaller values of /3 give an improvement 
of the fit over small sections of the light curve, but introduce sev- 
eral small spots that vary from one subset to the other Therefore, 
they are actually not required by the data. On the other hand, by 
increasing the entropy of the maps as much as possible, we ob- 
tain smoother maps with a regular variation of the spot pattern 
from one subset to the other, which indicates that we are prop- 
erly modelling the overall active region distribution on the star 
The two spikes of the best fit at HJD ~ 2454160 and ~ 2454168 
are due to the small variations of the spot configuration when we 
move from one data subset to the next. The reader may won- 
der about the possibility of improving the regularized best fit by 
considering smaller time intervals Aff. However, the improve- 
ment is negligible because the misfits at the matching points be- 
tween one subset and the next are a consequence of the signif- 
icant amount of regularization required by this light curve due 
to its relatively low information content. The effect is actually 
much larger in the present case than for CoRoT-Exo-2a, whose 
light curve had a sig nificantly greater in formation content than 
the present one (cf. iLanza et"ani2009l) . Therefore, decreasing 
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Fig. 2. Upper panel: The ME-regularized composite best fit to the out-of-transit light curve of CoRoT-Exo-4a obtained for Q - 4.5. 
The flux is in relative units, i.e., measured with respect to the maximum observed flux in the de-trended normal point light curve. 
Lower panel: The residuals from the composite best fit versus time. In both panels, the dashed vertical lines indicate the epochs of 
mid transits. 



Aff in the present case does not alleviate the problem, but pro- 
duces the opposite effect, owing to a greater number of matching 
points along the time series. 

The distributions of the spotted area vs. longitude are plotted 
in Fig. [3] for the seven mean epochs of our individual subsets 
adopting a rotation period of 9.20205 days. The longitude zero 
corresponds to the point intercepted on the stellar photosphere 
by the Hne of sight to the centre of the star at HJD 2454135.0918, 
i.e., the sub-observer point at the initial epoch. The longitude 
increases in the same sense as the stellar rotation and the orbital 
motion of the planet, which has a fixed subplanetary point on the 
star in the adopted reference frame. This choice allows a direct 
comparison of the active region longitudes with the subplanetary 
longitude, although one cannot correlate the active regions to the 
dips in the light curves in a straightforward way. 

Three main active regions can be identified in Fig. [3] and 
their migration has been traced with different straight lines. 
Specifically, their rotation rates with respect to the adopted refer- 
ence frame are found by a linear best fit, i.e., assuming a constant 
migration rate. They are listed in Table |2]for three values of the 
facular-to-spotted area ratio Q, because the longitudes of the ac- 
tive regions derived from our model have so me dependence on 
that parameter (cf. iLanza et"ani2007l l2009l) . The active region 
indicated as ARi corresponds to the one marked with a three- 
dot-dashed line in Fig. [3] AR2 to that marked with a long-dashed 
line, and AR3 to that marked with a dot-dashed line, respectively. 



Table 2. Relative migration rates AO/Q of the active longitudes 
for different facular-to-spotted area ratio Q and rotation period 
^"101 ■ 



Q 



(day) 



AR, 



ARo 



AR, 



0.0 9.202 

4.5 9.202 

4.5 8.870 

7.0 9.202 



0.106: 
0.108: 

0.063 : 

0.100: 



0.025 
0.008 
0.011 
0.011 



0.034 : 
0.052 : 

0.010: 

0.046 : 



0.008 
0.009 
0.010 
0.010 



0.070 ± 0.024 
0.100 ±0.024 
0.073 ± 0.036 
0.128 ±0.010 



The relative amplitude of the surface differential rotation, es- 
timated from the difference between the greatest and the lowest 
migration rates, is AQ/Q = 0.072 ± 0.026, 0.057 + 0.015, and 
0.0825 ± 0.014 for 2 = 0, 4.5, and 7.0, respectively. For the 
sake of completeness, we have also computed ME spot mod- 
els adopting Q - 4.5 and a rotation period of 8.870 days, 
that comes from the autocorrelation analysis of the light curve 
(Aigrain et al. 2008). Again, three main active longitudes are 
seen in the longitude distribution of the spotted area. Their mi- 
gration rates are listed in Table [2] and the derived relative am- 
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plitude for the differential rotation is AQ/O - 0.063 ± 0.038, 
in very good agreement with the resuh obtained for a rotation 
period of 9.20205 days. It is interesting to note that a rotation 
period of 8.870 days makes the migration rate of AR2, i.e., the 
active longitude with the greatest spotted area, virtually zero. 
This is in agreement with the fact that the autocorrelation of the 
light curve is maximized for that rotation period. 

Unfortunately, no information is available on the spot lati- 
tudes in our models, therefore, our AQ/f2 are only lower limits. 
These values suggest that CoRoT-Exo-4a has a surface differ- 
ential rotation comparable to that of the Sun, for which we can 
estimate a relative amplitude of 0.04-0.05 considering active re- 
gions confined within the sunspot belt, i.e., within + (35° - 40°) 
from the equator. 

The lifetimes of the active regions can be estimated from our 
spot models and range from 10 to 50 days, with longer lifetimes 
characterizing those having greater filling factors, as in the Sun. 

It is interesting to note that some active regions are located 
close to the subplanetary longitude during most of the time, al- 
though the highest concentrations of spotted (and facular) areas 
are around a longitude following the planet by about 100°- 120°. 
Conversely, the hemisphere centred around longitude » 240° has 
the minimum active region coverage (cf. Fig. [3]). The association 
between the subplanetary longitude and an enhanced spot activ- 
ity is particularly remarkable in the case of the models computed 
without facular contribution, i.e., for Q - (see Fig. |4|i. In six 
out of seven epochs, the ME models indicate a relative maximum 
of the spotted area within +50° from the subplanetary longitude. 
The time average of the spotted area in each longitude bin is 
plotted in Fig. |5]vs. longitude, together with the corresponding 
standard deviation. In the case with Q = Q (Fig.|5] lower panel), 
a well-defined relative maximum appears close to the subplane- 
tary longitude, with a height of ~ 3 standard deviations, which 
is unlikely to be due to a chance fluctuation. Moreover, the sub- 
planetary peak has a FWHM significantly smaller than that of 
the active longitude peak around 0° whose width reflects the mi- 
gration of its active regions along time. A greater migration rate 
of the active regions is also responsible for the smearing of the 
active longitude around 0° and the subplanetary longitude, when 
we consider the average of the models computed with Q = 4.5 
(Fig- El upper panel). 

The total spotted area varies only slightly versus time (see 
Fig. |6]) and may indicate a level of activity about 2-3 times 
that of the Sun at the maximum of the eleven-year cycle. It is 
important to notice that the total spotted area is determined here 
after removing the long-term decrease of the flux observed in the 
original light curve by Aigrainet al. (2008). Therefore, it does 
not show the long-term linear increase that one would expect in 
the case that the long-term trend were not removed. Moreover, 
the absolute value of the area depends on the adopted spot con- 
trast Cs and the value of Q. For instance, a lower spot temperature 
would imply a stronger contrast and thus a smaller area, but the 
relativ e variations of the area are largely independent of Cs and 
Q (cf. iLanza et alj|2009l) . 



6. Discussion 

Our results confirm that the rotation of CoRoT-Exo-4a is quasi- 
synchronized with the orbital motion of its planet. Minimum 
spot migration rates are obtained for a rotation period of 8.87 
days, which corresponds to the rotation period of the active lon- 
gitude with the largest spot concentration. However, even adopt- 
ing that period, the other active longitudes continue to display 
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Fig. 3. The distributions of the spotted area vs. longitude at the 
labelled times (HJD - 2450000.0) for Q = 4.5. The plots have 
been vertically shifted to show the migration of individual spots 
(relative maxima of the distributions) versus time. The vertical 
dashed lines mark longitudes 0° and 360°, beyond which the dis- 
tributions have been repeated to help following the spot migra- 
tion. The vertical dotted line marks the subplanetary longitude. 
The dot-dashed, three-dot-dashed and long-dashed lines trace 
the migration of the most conspicuous spots detected in the plots 
(see the text for details). 



a significant migration. Our interpretation for the different mi- 
gration rates of the active longitudes is that the corresponding 
active regions are located at different latitudes on a differen- 
tially rotating star Nevertheless, the different drift rates of the 
active longitudes may be produced by a phenomenon different 
from surface differential rotation. For instance, in the frame- 
work of a solar analogy, the rotation rate of sunspot groups 
during the early stages of their evolution (< 2 - 3 days from 
their appearance) is about 2 percent greater than that of recur- 
rent g roups and decreases steadily during their later evolution 
(e.g., IZappala & Zuccarellolll99l1 ). However, we believe that it 
is unlikely that a similar phenomenon may explain the results 
obtained for CoRoT-Exo-4a because our active longitude drift 
lasts for several tens of days. In conclusion, surface differential 
rotation in combination with active regions located at different 
latitudes seems to be the most likely explanation for the results 
presented in Sect.|5] Although the observed active longitudes ro- 
tate faster than the planet, the existence of a latitude exactly syn- 
chronized to the planetary orbit cannot be excluded because we 
have no information on the latitude range covered by the active 
regions. 

The quasi-synchronous rotation status of CoRoT-Exo-4a 
might have been attained when the star, although already close 
to the ZAMS, was still magnetically locked to the inner parts of 
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Fig. 4. The same as Fig.|3] but for Q = 0.0. 
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planet could not significantly change the angular moment um of 
the planet and the star, as noticed bv 'Aigrain e t alj (l2008f ). The 
magnetic braking of the star after the dissipation of the disc is 
expected to be negligible given it s mid-F spectral type and an 
initial rotation period of ~ 9 days jWolff & Simonll 19971) . 

The surface differential rotation found through our spot mod- 
elling is not remarkably dependent on the adopted value of the 
facular-to-spotted area ratio Q and is at least 6-7 times that 
derived for CoRoT-Ex o-2a that is rotatin g in 4.52 days and is 
of spectral type G7V (iLanza et alj|2009l) . For t B ootis, having 
Teff = 6360 K and a rotation period of ~ 3.3 davs. ICatala et al.l 
(12007) found a relative differential rotation AQ/Q « 0.15 - 0.18 
between the equator and the pole by means of line profile mod- 
elling. This compares well with our value, given that our deter- 
mination is only a lower limit. The same method was applied 
by Reiners (2006) to a sample of F-type stars finding projected 
pole-equator relative amplitudes that are comparable to that of 
CoRoT-Exo-4a when one considers objects with a similar ef- 
fective temperature and 15 < vsin/ < 30 km s Note that 
the small v sin ; of our star prevents the application of line pro- 
file analy sis methods to measure its surface differential rotation 
(lReinersll2006) . therefore, a spot modelling approach is the only 
possibility to estimate its surface shear. 

Theoretical predictions for the amplitude of the surface dif- 
ferential rotation in main-sequence F-type stars have recently 
been published by Kiiker & Riidiger (2005). Their model is in 
general agreement with the present findings although the lack of 
a precise theory for the interaction between stellar rotation and 
turbulent convection makes a detailed quantitative comparison 
premature. 

The possibility that stellar activity is influenced by the hot 
Jupiter in CoRoT-Exo-4a is suggested by the present models, in 
particular by those assuming that the photometric modulation is 
produced only by cool spots (i.e., with Q = 0). The plots in Fig.|5] 
suggest that the formation of spots may be triggered in an active 
region located at the subplanetary longitude, i.e., the emergence 
of magnetic flux may be promoted there in some way by the 
close planet. It is important to note that tidal effects cannot ac- 
count for this phenomenon because there is only one region of 
enhanced activity on the star during most of the time, and not 
two separated by 180°, as expected in that case. The possibil- 
ity that a hot Jupiter may affect in some way the dynamo action 
in the subsurface layers of the convection zone of its host star 
has been discussed by Lanza (2008 ), and the present results give 
some support to his speculations. However, a word of caution is 
certainly appropriate in the present case because the time interval 



Fig. 5. Upper panel: The spotted area per 18° longitude bin av- 
eraged over time versus longitude for Q - 4.5. The semiampli- 
tudes of the error bars are equal to one standard deviation of the 
spotted area in the corresponding bin. Lower panel: The same as 
in the upper panel, but for Q-Q. The vertical dotted line marks 
the subplanetary longitude, while the dashed vertical lines mark 
longitudes 0° and 360° beyond which the distributions are re- 
peated. 



a circumstellar disc. If the planet had already migrated to the ra- 
dius where the disc was truncated by the stellar magnetic field 
(Romanova & Lovelace 2006), the system started its evolution 
in a rotational and orbital configuration close to the present sta- 
tus. The subsequent tidal interaction between the star and the 




4150 4160 4170 4180 
Time (JD - 2450000) 

Fig. 6. The total spotted area as derived from the regularized ME 
models vs. time for Q - 4.5. 
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covered by our data is rather limited. Therefore, the supposedly 
enhanced activity related to the planet might still be the result of 
statistical fluctuations in the appearance of small active regions 
over the star Moreover, the relationship between enhanced ac- 
tivity and subplanetary longitude appears to be less clear when 
we consider models with Q > 0, also because the migration rates 
of the main active longitudes are greater. 

7. Conclusions 

Our analysis shows that the rotation of CoRoT-Exo-4a is quasi- 
synchronized, with the orbital motion of the planet. Specifically, 
the active longitude having the greatest spotted area shows an 
angular velocity that is only 3-5 percent greater than the or- 
bital angular velocity. Assuming that it traces the angular ve- 
locity of stellar rotation at some latitude, this translates into an 
upper limit for the difference between stellar and orbital angu- 
lar velocities at that latitude. Moreover, we find evidence for 
a significant differential rotation of the star, with a lower limit 
AQ/f2 = 0.057 ±0.015, when the optimal facular-to-spotted area 
ratio Q = 4.5 is adopted. 

It cannot be excluded that the planet may affect in some way 
the longitude distribution of the active regions on the photo- 
sphere of the star Specifically, there appears to be an enhace- 
ment of spot activity close to the subplanetary longitude, partic- 
ularly for the models computed without faculae. Unfortunately, 
the short duration of the present time series (only 58 days, i.e., 
6.3 stellar rotations) does not allow us to derive definite conclu- 
sions on such a kind of star-planet magnetic interaction. 

A long-term monitoring of the star from the ground may help 
to clarify this issue, e.g., through a measurement of the rotational 
modulation of the chromospheric Ca II H & K flux that should 
be a good indicator of a po ssible star-pla net interactio n, as in 
the ca se of t Boo or v And dWaUcer et al.li2008; .Shkolnik et al.1 
l2008h . 
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